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We present a technique for manipulating the nuclear spins and the emission polarization from a single 
optically active quantum dot. When the quantum dot is tunnel coupled to a Fermi sea, we have discovered 
a natural cycle in which an electron spin is repeatedly created with resonant optical excitation. The sponta- 
neous emission polarization and the nuclear spin polarization exhibit a bistability. For a a + pump, the emission 
switches from a + to o~ at a particular detuning of the laser. Simultaneously, the nuclear spin polarization 
switches from positive to negative. Away from the bistability, the nuclear spin polarization can be changed 
continuously from negative to positive, allowing precise control via the laser wavelength. 



PACS numbers: 73.21.La, 78.67.Hc 

Semiconductor quantum dots are very attractive for appli- 
cations as qubits [1] and sources of quantum light [2-4]. Ver- 
satile materials are the III-V semiconductors, notably GaAs 
which has established itself as the workhorse material. A sig- 
nificant property is that all the Ga, As, and In isotopes have 
large nuclear spins. In a typical quantum dot there is an inter- 
mediate number of atoms, too large to use each nuclear spin 
as a resource yet too small for efficient cancellation in the to- 
tal spin, and noise in the nuclear spins limits the electron spin 
coherence to just <~ 10 ns through the hyperfine interaction 
[5-7]. An emerging theme is that the nuclear spin noise may 
be reduced by narrowing the distribution [8-11] and that the 
nuclear spin ensemble may represent as much opportunity as 
trouble. Currently, schemes exist to tune both the optical tran- 
sition energy [ 1 2] and the selection rules [ 1 3] of a quantum dot 
in situ, but presently, the possibilities of using nuclear spins 
beneficially are limited. 

We present here a new control over the electron-nuclear 
spin interaction on driving an optical transition resonantly. 
Dynamic nuclear polarization at the single quantum dot level 
is established [14-18]. The crucial advance here is to operate 
in the tunneling regime [17, 18] where we discover a natural 
cycle. There are two interrelated features. First, spontaneous 
emission following resonant excitation either preserves the 
circular polarization of the source or inverts it. For instance, 
with a cr + pump, we can switch from predominantly a + to 
cr~ emission either with a small change in pump wavelength 
or device bias allowing the polarization of a single photon 
source to be controlled in situ. Secondly, the resonant exci- 
tation creates a large nuclear spin polarization which changes 
sign abruptly at the bistability, a new feature compared to the 
bistabilities following nonresonant optical excitation [14, 19]. 
At smaller laser wavelengths, the nuclear spin polarization 
changes monotonically from a large negative value to a large 
positive value. This bidirectional tuning is demonstrated here 
at low magnetic fields (0.5 T), and complements the optical 
dragging effect at high magnetic fields [12]. Control of the 
nuclear spins via the optical wavelength is a powerful route to 



narrowing the distribution [12] and to tuning the quantum dot 
exciton over tens of //eV. 

Our experiments use a field effect device in which InGaAs 
self-assembled quantum dots are in tunnel contact with an n + 
GaAs Fermi sea via a 25 nm thick GaAs tunnel barrier [20]. A 
voltage is applied to a Schottky contact on the sample surface, 
150 nm above the quantum dot layer, at 4.2 K. Photolumines- 
cence (PL) is excited either nonresonantly at 830 nm wave- 
length, or resonantly using 13 kW/cm 2 from a tunable narrow- 
band cw laser. The PL is dispersed with a monochromator and 
detected with a CCD array detector, a system with resolution 
50 /ieV. The polarization of excitation and collection are inde- 
pendently controlled. A small magnetic field, B z = +0.5 T, 
is applied along the growth z direction. 

Excited nonresonantly, the photoluminescence from a sin- 
gle quantum dot shows a clear charging step from the neutral 
exciton, X°, to the negatively charged trion, X 1 ^ . The en- 
ergies of the initial states \X°) and and their corre- 
sponding final states, |0) (vacuum) and |e) (single electron), 
as a function of gate voltage are shown in Fig. 1. In the fi- 
nal states (no hole present), the ground state charges from |0) 
to |e) at a more positive voltage than the change in the initial 
states from \X°) to \X X ~) [20], a consequence of the differ- 
ent electron-hole and electron-electron Coulomb energies. A 
"hybridization region" is created, a voltage region in which 
both excitons are tunnel coupled to the Fermi sea, X a in the 
initial state, X 1 ^ in the final state [20]. We show here that 
this region is ideal for controlling the electron-nuclear spin 
interaction. 

Figure 1 shows the result of pumping the |0) «-> \X°) tran- 
sition of a single quantum dot. Over a small region of volt- 
age, X x ~ PL is observed, redshifted by 6 meV with respect 
to the laser. A comparison with the nonresonantly excited PL 
demonstrates that this region corresponds to the low bias edge 
of the X x ~ plateau, i.e. the hybridization region, and that the 
resonantly excited PL has the X° energy. In terms of the level 
diagram in Fig. 1, the dot is initially in the vacuum state |0). 
The laser then creates an X°, which, although neutral, is un- 
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FIG. 1. Top: Photoluminescence (PL) at 4.2 K from a single quan- 
tum dot versus bias driven with excitation at the X energy. X 1 ^ PL 
appears in a narrow range of voltage, the hybridization regime. Bot- 
tom: Energy dependence versus bias for the quantum dot vacuum 
state |0) and the single electron state je), showing a crossing where 
the ground state changes. X° and X 1- cross at lower bias on ac- 
count of the hole. Within the hybridization region, automatic cycling 
takes place when a laser is tuned to the |0) <H> \X°) transition. An 
electron tunneling from the Fermi sea turns the \X°) into re- 
combination leaves the system in state j e) ; tunneling out returns the 
dot to |0). 



stable with respect to tunneling. Electron tunneling into the 
dot (time scale ~ 50 ps, considerably shorter than the radiative 
lifetime of ~ 1 ns) creates an X 1 " which then recombines. 
After spontaneous emission, the dot is in the |e) state. Now 
that the hole has disappeared, this state is also unstable with 
respect to tunneling: electron tunneling out of the dot (time 
scale ~ 10 ps) returns the dot to |0) whereupon the process 
can be repeated. This cycle offers a number of attractive fea- 
tures. First, the X° spin is determined by the polarization of 
the laser through the optical selection rules. Second, the cycle 
round-trip time is small, just ~ 1 ns, limited only by sponta- 
neous emission. Third, the redshift of the PL with respect to 
the excitation makes it easy to distinguish spontaneous emis- 
sion from scattered laser light even though one of the transi- 
tions is driven resonantly. The PL is useful in its own right 
as an antibunched source. It also provides an in situ monitor 
of the nuclear spin polarization through the Overhauser shift. 
Finally, the process can be described quantitatively with no ad 
hoc assumptions. 

The main experiment consists of monitoring the X 1 ^ PL as 
a function of laser detuning with respect to the X° transition 
for a constant pump polarization, e.g. a + . A PL spectrum 
is recorded for both a + and <r~ polarizations. These counts- 
energy spectra are fitted to Lorentzians [21], yielding both the 
signal energy E(a + /~ ) (center of Lorentzian), and the signal 
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FIG. 2. Left (right) panels: Experimentally measured (calculated) 
signal intensity, polarization degree, and Overhauser shift versus 
laser energy (laser detuning) for a a + pump and an external field 
of +0.5 T at fixed bias in the center of the hybridization region. In 
the experiment, the laser is tuned close to the |0) o \X°) transition 
and the dot (at 4.2 K) is the same as in Fig. 1. 



intensity S l (cr +//_ ) (area under Lorentzian). Figure 2 shows 
both S(cr + ) and S(a~) for a <r + pump, and the associated 
polarization degree 



P = 



S(a+)-S(a-) 
S(a+) + S(a~)' 



At large negative and positive detunings, the PL has largely 
a + character with P up to 0.76 ± 0.05. This is the intuitive 
result from the selection rules. Absorption of a a + photon 
with spin angular momentum +K along z creates an |f|\L) ex- 
citon consisting of a heavy hole f|~ with angular momentum z 
projection +f 7i and an electron J, with — An electron tun- 
nels in to form the X 1 ^ exciton, |f|"4t) ■ Hole spin relaxation is 
slow compared to recombination [22, 23] such that recombi- 
nation |f|"it) ~~ It) creates a a + photon. The counterintuitive 
result in Fig. 2 is that close to the center of the resonance, 
the PL has an inverted polarization degree, with P ~ —0.7. 
Strikingly, P changes abruptly at a particular detuning. 

An indicator that the nuclear spins are involved is provided 
by the Overhauser shift 



E{a+) - E(a-) - gxt&B, 
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A n is interpreted as an energy shift of the unpaired electron 
spin in the X 1 ^ final state arising from the nuclear spin po- 
larization along z. Its determination requires a knowledge of 
the exciton g-factor, and we measure gx = 1-55 ± 0.10 as 
described in [21]. Close to the center of the resonance we 
now find that A n switches sign exactly at the point where 
P switches sign. The Overhauser shift is related to the av- 
erage nuclear spin z projection (J z ) (in units of Ti) through 
A n ~ —A(I Z ) [21]. Taking the coupling constant A ss 
90 /J,eV, an averaged value for Ino.5Gao.5As [24], we find that 
(I z ) « +0.36 <H> —0.36. Full polarization corresponds to 
(I z ) = ±2.25, where / = 2.25 is the average nuclear spin 
quantum number in the dot. 

The abrupt jump in P corresponds to a bistability. With a + 
excitation, in state I (II) the dot emits a + (<j~) photons and the 
nuclear spins point up (down). The bistability is demonstrated 
clearly in the hysteresis curve of Fig. 3 (top). In this case, 
the laser energy was tuned in fine steps (less than 0.5 /ieV), 
blocking the laser path for about 30 s between each data point 
during which time the state of the system was always pre- 
served. At more positive laser detunings, the polarization de- 
gree P and the nuclear spin polarization are continuous mono- 
tonic functions of detuning, changing from large negative to 
large positive values. Correspondingly, A n goes smoothly 
from +30 to —35 ^eV. The total electron Zeeman splitting, 
g° fiQB z = g c /J,BB z + A(I Z ), changes sign at the bistability, 
followed by continuous tuning from —45 to +20 /ieV (tuning 
of effective electron g-factor g° s from —1.6 to +0.7). 

To switch from state I to state II, it is more convenient to 
change the gate voltage than the laser wavelength. We have 
achieved this by exploiting the Stark effect of the exciton en- 
ergy. Figure 3 (bottom) demonstrates controlled switching be- 
tween state I and II by applying voltage pulses to the gate, 
monitoring the state of the system via the a + PL. The system 
is initially in state I. It is forced into state II with a negative 
voltage pulse, equivalent to moving the laser energy up and 
back down again. This results in a lower a + PL, the signature 
of state II. Analogously, we can switch the system back into 
state I with a positive voltage pulse. In between these volt- 
age pulses, the laser is turned off. When it is turned back on 
again ~ 30 s later, the system always adopts its original state, 
demonstrating a slow I-II relaxation rate (< 0.1 s _1 ). 

We present a quantitative model to describe these results. 
The two crucial ingredients are first, a coherent coupling be- 
tween |f|- 1) and l-IJ-t). the so-called fine structure which arises 
from the anisotropic part of the electron-hole exchange, and 
second, a hyperfine coupling between the nuclear spins and 
the unpaired electron spin. A full description of the model is 
given in [21]. 

First, we calculate the effect of the laser field on the dynam- 
ics of a five-level system, consisting of the vacuum state |0), 
the two X° exciton states, |f|+) and |JJ-t), and the two X 1 ^ 
states, Ifl+t) an d 1-tJ-tD- The laser is <r + polarized and drives 
the |0) <-» Iffl) but not the |0) «-» |JJ-t) transition on account 
of the selection rules. The optical Rabi energy is hft, the de- 
tuning H6 — Tiui — Tlujq, where tu is the angular frequency of 
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FIG. 3. Top: Bistability between state I (a + PL: nuclear spins up) 
and state II (a~ PL: nuclear spins down). P and A n were measured 
as the laser was tuned. The laser was blocked for 30 s between each 
point. Bottom: Demonstration of switching with gate voltage pulses 
(40 mV, 5 s duration) by measuring the a + PL; strong PL signifying 
state I, weak PL signifying state II. In between voltage pulses, the 
laser was turned off. Both curves were recorded at +0.5 T, 4.2 K, 
using a a + pump and the same dot as in Figs. 1 and 2. 



the laser and huj is the eigenenergy of |f|+) and |JJ-t) without 
magnetic field. Coupling between |f|+) and |JJ.f) is character- 
ized by the fine structure fouj[ s . Decay processes are sketched 
in the level diagram, Fig. 4 (top). The neutral excitons can 
decay by spontaneous emission to |0) at rate Tq 1 ; or they can 
become trion states via tunneling at rate +T . Starting with 
the entire population in the ground state, we use the master 
equation for the density matrix to determine the occupation 
probabilities andp^^ of the trion states after time t\, 

the spontaneous recombination lifetime of X 1 ^, resulting in 
the rates of creating an f, i electron via optical recombination. 

After trion recombination, the free electron interacts with 
the N quantum dot nuclei through the contact hyperfine in- 
teraction before it tunnels out at rate r~^ t . The spin flip-flop 



probability ps is 



Pa 



27T 2 ut 
(47 + eK 2 ut + ft 



where 7 = ^ (/ - \(I Z )\) and £ = {g c ^B z + A{I Z ))\ 
with g c as the electron g-factor [21]. The combination of elec- 
tron creation rate and flip-flop probability results in a dynamic 
equation for the nuclear spin polarization. (I z ) is driven up 
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FIG. 4. Top: The five quantum states in the simulation showing an 
optical coupling (Rabi energy HQ,) and a coherent coupling (energy 
7kjf s ) between the two neutral exciton states. Decay processes are 
drawn with dashed lines. Bottom: The calculated nuclear spin dy- 
namics as a function of laser detuning for a er + pump, +0.5 T exter- 
nal field, and the parameters described in the text. The solid (dashed) 
line shows the stable (unstable) solution for 4z{Iz) = 0. 



in the |^t) state, leading to electron spin I creation following 
tunneling in and recombination. When the laser is then tuned 
further, the allowed |0) o transition takes over and the 
cycle results in the creation of electron spin f. The creation 
of a particular electron spin leads to nuclear spin polarization 
which alters the energies of the |-l|t)> 1114) states via the Over- 
hauser field. This feedback results in a bistability close to the 
forbidden transition and continuous tuning thereafter. 

We have explored some of the parameter space theoreti- 
cally. For parameters close to the ones used in this experi- 
ment, a region of bistability exists when Fi ca k is small enough. 
A bistability is definitely possible even at zero magnetic field, 
provided that ri Ga k < 1 s _1 and that the tunneling times are 
increased relative to those in this experiment. The inversion in 
polarization can be enhanced to at least P = +0.85 O —0.85, 
again by increasing the tunneling times and also by optimizing 
the Wf s : f2 ratio. Furthermore, at detunings larger than those 
at the bistability, these parameters allow continuous control 
of the l-U-t), |1T4) eigenenergies from +40 to ±50 /xeV, and, 
following Refs. [21, 25], a reduction in the variance of the nu- 
clear spin distribution by factors ~ 5. All these features are 
attractive for spin qubits and single photon emitters. 

We thank Bill Coish for very helpful discussions and 
acknowledge funding from Swiss NSF, DARPA QuEST, 
EPSRC, The Royal Society (BDG), and ANR QUAMOS 
(BU). 



depending on p^^-), down depending on P\^i), and decays 
in the absence of driving with rate ri ea k: 

j^ 1 *) - j^r [p\u\) -P\m)]t = n ~ r icak(^ 2 )- 

We solve this equation numerically to find stable values of 
(I z ) as a function of laser detuning fid. At each solution one 
can also calculate the Overhauser shift A n and the polariza- 
tion degree P in the quantum dot emission [21]. 

Parameters are set by in situ characterization and by com- 
parison with previous experiments, making small tweaks to fit 
the experimental data in Fig. 2. We use the following val- 
ues [21]: TiVt = 23 /xeV, fru>[ s = 40 /ieV, To = 0.75 ns, 
n = 0.95 ns, r in = 35 ps, r out = 5 ps, N = 8.5 x 10 4 , 
ricak = 0.1 s -1 and g e = —0.5. Fig. 4 (bottom) contains a 
plot of 4j.(I z )< showing that the solution for (I z ) changes from 
positive to negative with a region of bistability. The calculated 
P and A n are plotted in the right panels to Fig. 2. Close to the 
optical resonance, there is an excellent agreement with the ex- 
perimental results. 

The theory offers an explanation for the counterintuitive in- 
version of the PL polarization. When the er + -polarized laser 
comes into resonance with the forbidden |0) O |-IJ-t) transi- 
tion, a combination of the allowed |0) ■<-» transition and 
the Iff!) o \iyi) coupling causes the population to build up 
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In "Controlling the Interaction of Electron and Nuclear Spins in a Tunnel -Coupled Quantum Dot" [1] we report a technique 
for manipulating the nuclear spins and the emission polarization from a single optically active quantum dot. In order to describe 
this system theoretically, we divide the calculation into two interlinked parts. Process I: the laser field induces dynamics in a 
five-level electronic system. The final event is X 1 ^ recombination. Process II: nuclear spin polarization builds up via electron- 
nuclear spin flip-flop processes. Process I provides an unpaired electron spin as input to process II. Conversely, process II creates 
an Overhauser field which is an input to process I. In the presence of nuclear spin leakage, a master equation for the nuclear spin 
polarization (I z ) has stable solutions where d(I z )/dt = 0. Associated with the stable solutions for (I z ) are particular values of 
the Overhauser field, and hence particular values of the photoluminescence energies and intensities for the two light polarizations, 
(7+ and <t~. The details of this calculation are described here along with the input parameters. The results of the calculation can 
be compared directly with the experimental results (Fig. 2 of [1]). Furthermore, we comment on the consequences of the theory 
for the fluctuations of the nuclear spin ensemble. 



I. COHERENT EVOLUTION OF THE FIVE-LEVEL SYSTEM 



The dynamics of the five-level quantum system induced by the laser field are calculated with the density matrix. The ground 
state is the empty dot, |0). There are two neutral exciton states, |f|\|,) (hole spin up, electron spin down) and (hole spin 
down, electron spin up). Defining the ground state energy as zero, the diagonal terms of the Hamiltonian matrix H are 

(U\ H \U) = Kw + 1/2 • ( 9 xVbB z - A(I Z )) , 
(n\ H \n) = Klj - 1/2 • (g x p B B z - A{Q) , 

where Uojo is the eigenenergy of \^\.) and |J|t) m the absence of a magnetic field. Coupling to the external magnetic field B z 
is determined by the exciton g-factor gx- The electron also interacts with the internal (Overhauser) magnetic field, given by 
the expectation value of the nuclear spin z projection (I z ) (in units of K) and the averaged effective coupling constant A. The 
<7 + -polarized laser of angular frequency uj = lu + 5, with Ti5 as detuning, drives the 1 0) <H- |f|-|) but not the |0) j-IJ-t) transition 
on account of the selection rules. As the laser is always close to the resonance, the rotating wave approximation applies. The 
off-diagonal elements in the Hamiltonian are 

(0\H\U) =hn/2-e^, 
in\ H \U) = W2, 

where TiVl is the optical Rabi energy, and Tioj{ s is the fine structure splitting arising from the anisotropic part of the electron-hole 
exchange [2]. The two exciton states can decay either by spontaneous emission to the ground state at rate Tq , or by tunneling 
of an electron into the quantum dot at rate T;" 1 . Tunneling converts \^\\) into the X 1 " state |f|-|t); l-H-t) mt0 tne other X 1- state 
l-llti)- The trion states decay by spontaneous recombination, at rate r-f 1 , to states \\), |4-), respectively. Figure 4 (top) of [1] 
summarizes all the coherent couplings and decay processes. 

At t = 0, the entire population resides in the ground state. The probability of finding the system in state |s) at times t > 
corresponds to element (s\ p \s) of the density matrix p, and we use the master equation with the decay processes in the Lindblad 
form, 



d „ i 

dt p = ~h 



H,p 



+ Dp, 



Dp = Y^ r im>^i«> (\ n ) H p \ m ) ( n \ - \ \ m ) H p - \p \ m ) H j i 

m,n ^ ^ 

to calculate its coherent evolution [3]. The dissipator Dp accounts for the decay processes, where F| m ^|„^ is the transition rate 
from state |m) to |n). 
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Of particular interest are the trion populations, p^-f) = (itit\ P I It-It)' P|4ti> = (4J-T4-I P l-D-t4-)> since they determine both the 
X 1 ^ emission intensity and the probability of creating a conduction level electron with spin up (down) after trion recombination. 
The spontaneous X 1 ^ recombination time is t\, and we hence take {p\^^) /Ti)\ t=Tl and (p\^\) / T-\)\ t=Tl as the rates of creating 
an f, 4- electron, respectively. This approximation, essentially that the dynamics approach the steady state faster than spontaneous 
recombination, is well justified for laser energies near resonance. For large detunings, the approximation works less well, but 
going beyond this approximation would greatly complicate the calculation at the expense of transparency. Figure 1 (left) shows 
the evolution of the trion population as calculated in the five-level system for a laser energy in the region of the bistability. The 
parameters are the same as used in the main article [1], summarized in detail in section V below. For positive (negative) nuclear 
spin polarization, the |ltit) (l-IJ-tt)) state i s preferentially occupied, leading to a + (a~) emission and dynamic nuclear spin 
polarization in positive (negative) z direction. 
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FIG. 1. Left: Simulation of the trion occupation probabilities pi^-f) and pi^fj.) in the five-level system of ground state, neutral and negatively 
charged excitons, taking the entire population in the ground state at time t = 0. The parameters are the same as used in the main article [1], 
summarized in detail in section V, with the detuning hS fixed in the region of bistability. Depending on the value of (I z ), either |lHt) or I-D-T4-) 
is preferentially occupied after the X 1 ^ recombination time Ty, determining the rate of creating an f, J, electron via optical recombination. 
Right: Measured counts-energy spectra with the laser tuned to the region of bistability. The dot, at 4.2 K, is the same as in Figs. 1-3 of [1] and 
the energy of the detected signal corresponds to X 1- . The experimental data are fitted to Lorentzians, yielding both the energies E(a + /~) 
(center of curve) and the signal intensities S(a + ^~ ) (area under curve). A clear shift in the energies is observed, arising from different nuclear 
spin polarizations (I z )- As expected from the theory (left), the a + signal is pronounced when the nuclear spin polarization is positive, and 
vice versa. Both theory and experiment use a a + pump and an external field of B z = +0.5 T. 



II. FLIP-FLOP OF ELECTRON AND NUCLEAR SPIN 

At time t = t±, after X 1 ^ recombination, the quantum dot is in the |e) state and contains a free electron. This electron 
interacts with the N nuclear spins in the quantum dot before it tunnels into the Fermi sea at rate T~^ t . Since the conduction level 
electron has an s-type Bloch function, its coupling to the nuclear spins is well described by the contact hyperfine interaction with 
effective Hamiltonian 

H h f = vaJ2 A k\^k)\ 2 Ik-S. 

k 

Here A 7 are the effective hyperfine coupling constants for atoms of type j, v denotes the volume per atom, ip(Rk) is the electron 
wave function at nucleus k, and /, S are the nuclear and electron spins in units of h [4]. We assume homogeneous coupling, 
which yields 

a result used in the diagonal elements of the five-level system. Based on Ino.5Gao.5As, we take A ~ 90 peV as the averaged 
effective coupling constant [4] and / = 0.75 • 3/2 + 0.25 • 9/2 = 2.25 as averaged nuclear spin quantum number. Adding the 
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Zeeman term H cxt , with g c as electron g-factor, the total Hamiltonian reads 

ffcxt + Hhf - 9*HbB x S x + A(I Z )S Z + ( T k S ~ + ^^ + ) . (2) 

k 

where 7^ = 1% ± i/£ and S ± = ^ ± zSj, are ladder operators, such that the final term describes a flip-flop of electron and 
nuclear spin. In the following it is assumed that the electron has spin up; the calculation for spin down is analogous. 

To estimate the flip-flop probability as a function of net nuclear polarization, we calculate the time evolution in the effective 
2D system {|M,t), \M + 1,4-)}, where \M) represents the ensemble of nuclear spin states with z projection M = N{I Z ). 
With the system in state \M, t) at time t' = t — t\ = 0, one finds from the von Neumann equation [3] that the population of 
\M + 1,1} oscillates in time according to 



P|M + u) = 4 ^sm {-^—t), (3) 



where 



(M + l,l\H h{ |M,t)| 2 , 
£ = | (M + 1,4.| H cxt + H M \M + 14) - (M, t| H cxt + H M \M, t) |* . 

From Eq. (2), £ = (g c ^B z + A(I Z ) ) 2 . Simple results for 7 only exist when the nuclear spins are either fully correlated or fully 
uncorrected, neither of which is the case here. With 7 = (I — we take an average of the two extremes, as derived 

in the appendix. Integrating Eq. (3) over the survival distribution, characterized by the tunneling time T out , finally yields the 
flip-flop probability 



f ^ e-'/^ 2 7 r Q 2 ut 

P« = / dt P\M+l.i) = — T2- 

J Tout (4 7 + T D 2 ut + K 2 



t'=0 



Referring to previous work on dynamical nuclear polarization, the result is of a rather standard form, with r ou t as correlation 
time [5, 6]. Since we operate at rather low magnetic fields, where the "bright" |f|-|) , |JJ-t) X° states are split from the "dark" 
llTt) j I-H4) X° states by far more than a hundred /zeV [2, 7], spin flip-flops among electron and nuclei are negligible in the 
presence of a hole. 



III. NUCLEAR SPIN DYNAMICS AND STABLE SOLUTIONS 



The combination of the electron spin creation rate and the flip-flop probability results in a dynamic equation for the nuclear 
spin polarization (I z ). It is driven up depending on pi-fj-m, down depending on and decays in the absence of driving with 

rate ri oa k. With ri(<~ ns) 3> T out (^ ps), the cycle round-trip time is simply n and the dynamic equation reads 

J t Vz) - [P\Ut) -P\m)] t = n - r icak(/,)- (4) 

The flip-flop probability pg is a function of (I z ), the populations P\^), P|jj.ti> are functions of both nuclear spin polarization 
and laser detuning. We solve Eq. (4) numerically to find stable values of (I z ) as a function of detuning M. Figure 2 plots the 
dynamic equation (4) as a function of nuclear spin polarization, parameters as in section V, as the laser is tuned through the 
bistability. Stable nuclear spin polarizations are found where ^j(I z ) crosses zero with negative slope. The stable nuclear spin 
polarizations are directly related to the Overhauser shift A n observed in the experiment. Comparing the X x ~ emission energies 
E(ct + /~) measured for cr + /~ polarization, 

E(a+) - E(a-) = g^ B B z + A n = gf^B z , 

where we take A n = —A(I Z ) according to Eq. (1). Figure 1 (right) shows four counts-energy spectra measured in the region of 
bistability and fitted to Lorentzian curves, where E(a + /~) are given by the centers of the respective fit functions. It illustrates 
a key feature of the presented scheme, the tuning of resonance energies via control over the net polarization in the nuclear 
spin bath. In [1], a change from —45 /ieV to +20 /ieV in the total electron Zeeman splitting g° s /j,rB z — g c [i^,B z + A(I Z ) is 
observed, so that the sign of the effective electron g-factor g° ff can be inverted and, using the continuous change in (I z ) at more 
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FIG. 2. Simulation of 4r{I z ), Eq. (4), with the parameters as summarized in section V and B z = +0.5 T. Stable nuclear spin polarizations 
are found where the horizontal axis is crossed with negative slope. For hSi = —42 /ieV (green), one solution exists at positive (I z )- For 
T1S2 = —38.5 /ieV (red), an additional stable solution exists at negative nuclear spin polarization, leading to a bistability. The solution at 
positive (I z ) disappears when the laser detuning is further increased, h$3 — —36 /ieV (blue). The theoretical width of the bistable region is 
2 ^ieV, and in general depends on the various input parameters. Widths > 10 fieW are achievable. 



positive laser energies, tuned precisely to zero. Accordingly, the effective exciton g-factor is tuned from 2.6 to 0.3, and 
tuning <?|? through zero should easily be possible with an optimized cycle. 

The Lorentzians fitted to the experimental data not only yield the line energies (center), but also the line intensities (area under 
the curve), denoted as S(<r + /~). In terms of the theoretical model, the <j + (cr~) intensity corresponds to the rate of creating 
an j" (\) electron via X 1 ^ recombination. For each stable solution of 4i{Iz) = we can hence calculate the er + /~ emission 
intensity and the associated degree of polarization, 



S{<7+) -S(a ) p\ m) -P\ m) 



S{(T+) + S(a ) p mt) +p lm) 



t = n 



IV. CALCULATION OF NUCLEAR SPIN DISTRIBUTION WIDTH 



The model also provides insight into the width of the nuclear spin distribution. The typical fluctuating nuclear magnetic 
field seen by the electron spin via the hyperfine interaction is ~ A/ {g c ^B y/N) [8, 9]. This corresponds to a fluctuation of the 
nuclear spin polarization (I z ) by typically 1/yN and, in terms of the associated Gaussian distribution function, a variance of 
<j 2 ~ 1/N. To calculate the nuclear spin distribution in the presence of the laser field, we follow the steps as summarized in 
Ref. [10] (supplementary information). Starting from a basic rate equation, reformulation in the continuum limit results in a 
Fokker-Planck equation which is solved for stable nuclear spin distributions. These solutions are of the form e F ^ I ~^\ expanding 
the argument F to second order in (I z ) around the stable nuclear spin polarization I z ,o yields the variance of the distribution. 
One finds 

where T t is the total spin flip rate divided by N, 

r * = [p\uv +mn)]t = ti + r icak|(4)l + 2r loak . 

The additional 2ri ca k accounts for intrinsic diffusion and ensures that a 2 = 1/N in the absence of a laser field. The simulation 
shows that the cycle in the particular experiment in [1] results in a reduction of the variance by a small factor (~ 1.7). However, 
the scheme is capable of narrowing the variance of the nuclear spin distribution by factors ~ 5 (similar to the numbers reported 
e.g. in [1 1]) by increasing the tunneling times slightly. 
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V. INPUT PARAMETERS 



The parameters for the simulation are set by in situ characterization of the quantum dot, comparison with previous experiments 
on the same sample, and by making small tweaks to fit the experimental data presented in Fig. 2 of [1]. The main features - an 
abrupt inversion in the polarization of both photoluminescence (PL) and nuclear spins, followed by continuous tuning through 
zero - are not sensitive to the parameters in the calculation. Reproducing the exact experimental results is however sensitive to 
the exact parameters, in particular fiujf s , T m , r ou t. an d Tioak- 

• The Rabi energy fofl is proportional to the oscillating electric field and can therefore be estimated from the power density 
at the quantum dot. This gives us a rough estimate of Ml A more precise determination is taken from [12] in which the 
Rabi energy is measured directly from the Autler-Townes splitting in a pump-probe experiment. This experiment used the 
same microscope and a sample from the same wafer as the present experiment. Given a laser power of 13 kW/cm 2 and 
the values reported in [12], we determine TiVl — 23 /ieV. 

• The fine structure splitting Tiujf s is highly scattered from dot to dot in this sample [13], 40 yueV is taken here, largely by 
fitting the results of the model to the experimental data. 

• For the exciton recombination times r (X°) and n (A" 1_ ), we take r = 0.75 ns and n = 0.95 ns from lifetime 
measurements of quantum dots in this wafer [14]. There is a systematic increase in the recombination lifetimes on going 
from X° to X 1 ^, and a systematic dependence on the X° exciton energy. Both these factors are included. 

• The time T; n = 35 ps for tunneling into the dot is estimated from r and the ratio of A" , X x ~ PL intensities in the 
hybridization region. The tunneling time for tunneling out of the dot is smaller as the electron is now above rather than 
below the Fermi energy (where the energy barrier is thinner). Experiments probing the spin cotunneling rate [15] and the 
photoluminescence in the hybridization region [16] on dots similar to the one used here point to a r out time of <~ 10 ps. 
Fine tuning by fitting the experimental data gives r out = 5 ps. 

• N, the number of quantum dot nuclear spins, is estimated from the extent of the electron wave function. The ground state 
for harmonic confinement, 



ipi(r) = Cie 



1 I x 2 I y I z 2 

2 \P X 11 ll 



is replaced by a step function ip% of amplitude C 2 = C\/^/e, the value where \ipi\ 2 4irr 2 dr, with r = y 'x 2 + y 2 + z 2 , 
peaks in the simple case of spherical symmetry. The lengths L- L are related to the quantization energies fkVi and effective 
electron mass m* via l t — y/Ti/ (m*w,). With hcj Xty — 30 meV [17] and hiu z — 95 meV [18] we calculate N — 8.5 x 10 4 
as the number of nuclear spins covered by ip 2 - This result is consistent with a measurement of the fluctuations in the 
Overhauser field [19] (see Section IV). 

The electron g-factor is g e = —0.5 as determined by optically detected single electron spin resonance performed on a dot 
in a very similar sample to the one used here [19]. For the exciton, gx = 1-55 was measured by assuming that at large 
detunings, where E(a + ) — E(a~) takes on a constant value, A n changes sign but not magnitude on switching the pump 
polarization, 

9X " 2^B Z (t E{a+) E{U ~^ + Pump + - E ^a- pump 



• The order of magnitude of the spin depolarization rate ri oa k results from direct experimental observations. As demon- 
strated in [1], the system always preserves its state when the laser, tuned to the region of bistability, is turned off and back 
on after times ~ 30 s. For times > 1 min this is no longer true, pointing towards leakage rates of <~ 0.01 to 0.1 s _1 . 
This measures the decay of the nuclear spin polarization with the dot in the vacuum state, |0). Optical pumping in the 
hybridization region causes the dot, averaged over time, to be partially occupied such that the relevant value of ri ca k is 
likely to be larger. We take r loak = 0.1 s _1 here as it gives the best fit to the experimental data. 
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Appendix: Flip-flop term in an effective 2D system 

Since (I z ) is observed via the Overhauser shift A n = —A(I Z ), we calculate the flip-flop probability as a function of nuclear 
spin polarization. Letting the total nuclear spin z projection, in units of h, be M, 

M = ^r k = N(i z ), 

k 

then the state of the nuclear spin bath is a linear superposition of states \J, M, \ J M ), where M is fixed, \M\ < J < NI is the 
quantum number for the total spin length, and \ J M is an additional quantum number whose range depends on M and J, such that 

NI J 

E E Ei j ' m '' a m'X j ' m '' a m<i = i- 

J=0(i)M'=-J A 

Assuming the free electron has spin up, the Hamiltonian _ff oxt + i?hf given in Eq. (2) couples the states \J, M, \ J M , t) with states 
| J, M + 1, Aj^ +1 ,4-), where the electron spin was simply added to the notation. For given M, considering the number of ~ 10 5 
nuclear spins, the system of interest is therefore high dimensional. However, it can be reduced to an effective 2D system. 

Since the energy depends only on the electron spin and the total nuclear spin z polarization, all states associated with |t), \i) 
have the same eigenenergy E\, E 2 respectively. The population at time t' = is distributed among the states \J, M, A^,t), 
which we formally replace by a state |M,t) with 100% occupation probability and eigenenergy E\. This state |M,t) is 
coupled to the states with nuclear spin z projection (M + 1) via the flip-flop term in the contact hyperfine Hamiltonian, 
Eqs. (1, 2). For the estimate of the spin flip probability, it only matters whether the system occupies state \M, f) or any of 
the other states | J, M + 1, X'[ I+1 ,l). Therefore the dynamics can equivalently be calculated in a reduced system, where the 
I J, M + 1, \ J M+1 , 1) are replaced by the effective state \M + of eigenenergy E 2 and effective coupling 7 where 



(M+ui£hf|M,t) = E E|< J ' M+1 ' A M+i'^f wa) 

J=\M+1\ A 



(A.l) 



The value of 7 strongly depends on the system. In the following, the two extremes of a strongly correlated and a fully uncorre- 
lated system will be considered. 



Strong correlation 



As a first limit, we assume a highly coherent system free of further interactions, such that transitions J «-> J+ 1 are impossible. 
The relevant flip-flop term in the hyperfine Hamiltonian is 




where J + = J2k^k acts as a ra ising operator for the total spin. Its effect on a state | J, M, T) (A irrelevant) within \M, t) is 

Aj+5-| J)M)t ) = A^(j _ M )( J + W+T) \J, M (A.2) 

which breaks down as soon as the expression under the square root is zero. Starting from an unpolarized system, the typical 
order of magnitude for the starting value of J can be estimated from the infinite temperature limit: 

j(j + 1) = (j- j) = E ^ • + E^ • ^> = E^ • = Ni v + !)> 

k^k' k k 

J~VNy/l(I + l). 

Given this value, Eq. (A.2) typically approaches zero when |(7 Z )| ~ + 1)/N. For N = 8.5 x 10 4 and I = 2.25, these 

so-called "dark" states would therefore be reached at nuclear spin polarizations of less than 1% in both positive and negative z 
directions. This is in clear contradiction with the experiment, where polarizations of <~ 20% were measured and dark states were 
not observed. 
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No correlation 



In the other extreme, we take an uncorrelated system without transverse spin coherence, so that there is no correlation among 
the x, y components of the nuclear spins. Starting from Eq. (A.l), 



NI 

T= E E|< J > M + 1 ' A M+i>+l^hf|M,t> 

J=\M+l\ A 

NI 



A? 
AN 2 



NI 



J2 J2\{J,M + l,X J M+1 \J2l+\M) 



J=\M+l\ A 



A 2 
AN 2 



Ew* E Ei j ' m + 1 ' a m+i)( j ' m+1 ' a m+iI ^\ M ) 



J=\M+l\ A 



AN 2 



k^k' 



equivalent to 1 , only non-zero contribution 
A 2 

AN 2 



1 ^<M|/ fc -J+|M) + -^]T(M|J fe -/+|M> 



A 2 
AN 2 



A 2 



J2 (M| (I 2 + / - {I z k f - II) \M) ^ —(I 2 + I- (I z ) 2 (Q) . 



This calculation represents the nuclear spin system when coherence in the nuclear spin bath is lost between two cycles via various 
mechanisms such as the dipole-dipole interaction among the nuclear spins. However, these interactions are weak and the cycle 
round-trip time is very short, <~ 1 ns, so that this limit isn't realistic either. 

Considering I = 2.25, we therefore take 7 = ^ (7 — | (I z ) |) as a compromise of the two extremes. An analogous calculation 
can be carried out when the electron spin is initially down. It is clear that 7 must decrease as | (I z ) | increases. The exact functional 
form of 7 is not crucial to describe theoretically the main experimental phenomena but it is important to describe the experimental 
results quantitatively. 
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